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LIS: Entry 2 of 3 



File: USPT 



Dec 10, 2002 



DOCUMENT -IDENTIFIER : US 6493632 Bl 

TITLE: Water saturation and sand fraction determination from borehole resi s t i v i t y 
imaging tool, transverse j_ nc juction logging and a tensorial water saturation model 



Brief Summary Text (15) : 

U.S. Pat. No. 5,656,930 issued to Hagiwara discloses a method of determining the 
h9.jrA.£?.5£.^ and tne anisotropy coefficient of 

a subterranean formation by means of an induction type logging tool positioned in a 
deviated borehole within the subterranean formation. In a preferred implementation, 
the induction type logging tool is first calibrated to determine a proportionality 
constant. A predetermined relationship between the proportionality constant, the 

phase shift derived rjesistivity_,__ the attenuation derived .resistivity^ the 

horizontal resistivity, the vertical resistivity, and the anisotropy coefficient is 
then generated and stored in the memory of a programmed central processing unit. 
During an induction logging operation, the phase shift derived resistivity and 
attenuation derived resistivi ty are then received and processed by the programmed 
central processing unit in accordance with the predetermined relationship to 
generate the horizontal resistivity, the vertical resistivity, and the anisotropy 
coefficient. These measured values of horizontal and vertical resistivities when 
combined with a predetermined relationship between the horizontal resistivity, the 
vertical res i s t i y 1 1 y , the net/gross ratio, and the ratio of the sand layer 
resistivity to the shale layer_res^ti_vity make it possible to obtain a net/gross 
ratio. However, there are many laminated reservoirs in which the sands may include 
dispersed shales. Interpretation of formation water saturation in such reservoirs 
can be in error if the combined effects of laminations, dispersed shales within the 
sand, and possible intrinsic anisotropy of the shales is not considered. 

Drawing Description Text (7) : 

FIG. 5 gives the steps of one subprocess of an embodiment of the invention for 
determination of water saturation from measured values of vertical an horizontal 
resistivi ties . 

Detailed Description Text (10) : 

With flat dips, the sinusoids have essentially zero amplitude. In one aspect of the 
present invention, the resistivity measurements are averaged circumf erentially and 
vertically within each identified layer to give an average resi stivi ty measurement 
for each layer identified above. Once this is done, the subsurface may be 
characterized by a number of plane layers, each of which has a constant 
resistivity. With the resolution of the button-electrode tool, these layers may 
range in thickness from a few millimeters to a few centimeters. 

Detailed Description Text (13) : 

As would be known to those versed in the art, a finely laminated sequence of layers 
having different resistivities exhibits a transverse isotropy on a larger scale 
where the wavelength of the electromagnetic wave is much greater than the layer 
thickness. This condition is easily satisfied even for propagation resistivity 
tools that, e.g., operate at a frequency of 2 MHz (with a 

wavelength . lambda .. apprxeq . 6 meters); for induction logging tools that have 
frequencies of the order of 50 kHz to 200 kHz, the wavelengths are even longer. For 
such interrogating frequencies, the layered medium is characterized by a horizontal 
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resistivity R.sub.h * and a v e r t; i c a 1 r e s i s t i v i t y R.sub.v * given by: ##EQU3## 
Detailed Description Text (17) : 

Referring again to FIG. 5, an induction or wave propagation tool is used to make 
measurements of the vertical and horizontal resistivity of the earth formations 
104. For example, U.S. Pat. No. 5,781,436 to Forgang et al , the contents of which 
are fully incorporated here by reference, discloses a method an apparatus for 
making measurements of horizontal and vertical resistivities of a transversely 
isotropic formation. 

Detailed Description Text (18) : 

The method disclosed by Forgang et al comprises selectively passing an alternating 
current through transmitter coils inserted into the wellbore. Each of the 
transmitter coils has a magnetic moment direction different from the magnetic 
moment direction of the other ones of the transmitter coils. The alternating 
current includes a first and a second frequency. The amplitude at the first 
frequency has a predetermined relationship to the amplitude at the second 
frequency. The relationship corresponds to the first and the second frequencies. 
The method includes selectively receiving voltages induced in a receiver coil 
having a sensitive direction substantially parallel to the axis of the 
corresponding transmitter coil through which the alternating current is passed. A 
difference in magnitudes between a component of the received voltage at the first 
frequency and a component of the voltage at the second frequency is measured, and 
conductivity is calculated from the difference in magnitudes of the components of 
the received voltage at the two frequencies. The Forgang patent is cited only by 
way of example of an induction device for obtaining horizontal and vertical 
r^istivji^ties of a formation and there are other teachings on obtaining these 
properties of subterranean formation. 

Detailed Description Text (19) : 

An example of a propagation resistivity tool for making measurements of horizontal 
and yertica is described by Rosthal (U.S. Pat. No. 5,329,448) 

discloses a method for determining the horizontal and vertical conductivities from 
a propagation logging device. The method assumes that .theta., the angle between 
the borehole axis and the normal to the bedding plane, is known. Conductivity 
estimates are obtained by two methods. The first method measures the attenuation of 
the amplitude of the received signal between two receivers and derives a first 
estimate of conductivity from this attenuation. The second method measures the 
phase difference between the received signals at two receivers and derives a second 
estimate of conductivity from this phase shift. Two estimates are used to give the 
starting estimate of a conductivity model and based on this model, an attenuation 
and a phase shift for the two receivers are calculated. An iterative scheme is then 
used to update the initial conductivity model until a good match is obtained 
between the model output and the actual measured attenuation and phase shift. 

Detailed Description Text (20) : 

The horizontal and vertical ....resi^^iy^ties^ obtained at 104 are then inverted to give 
a layered model o f r e s i a it iyities 106. U.S. Pat. No. 5, 854, 991 issued to Gupta et 
al, the contents of which are fully incorporated here by reference, discloses such 
a method for inversion of transverse electromagnetic induction well logging 
measurements. Co-pending U.S. patent application Ser. No: 09/052,462 having the 
same assignee as the present application, and the contents of which are fully 
incorporated herein by reference, discloses a method of inversion of data from 
propagation resistivity tools. The f 462 application also discloses an ambiguity in 
the inversion of the resistivity measurements and a method for dealing with the 
ambiguity. 

Detailed Description Text (21) : 

In the inversion method used by Gupta et al, model is generated of the axial 
distribution of the horizontal and vertical conductivities, from induction signals 
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acquired by the instrument using two-frequency alternating current. The model is 
generated by calculating an initial estimate of the conductivity distribution and 
axially inverting the estimate with respect to the measurements made using the two- 
frequency alternating current. Shoulder correction is applied to measurements made 
by the instrument using single-frequency alternating current. An estimate of the 
radial distribution of the conductivities is generated from the shoulder corrected 
induction signals acquired using the single- frequency alternating current. A 2- 
dimensional model is made of the conductivity distribution from the model of axial 
distribution and from the estimate of radial distribution. In one embodiment of the 
present invention, the initial model for the inversion is based at least in part on 
data acquired by the resistivity imaging tool 105. In particular, while the 
£?.sis.t:.iYAty. imaging tool may need to be normalized in some way to correct the 
resistivity measurements, the layer boundaries determined by the resistivity 
imaging tool serve as a good starting point for the layers used in the inversion of 
the transverse induction logging tool data. In another embodiment of the invention, 
the initial layers for the model may be determined from other high resolution 
logging tools, such as a Laterolog . TM. or a gamma ray logging too. 

Detailed Description Text (23) : 

Once the inversion is performed, the horizontal and vertical resistivities obtained 
therefrom are analyzed using a petrophysical model 108. In one embodiment of the 
invention, the values of vertica_l__a_^ thus obtained are 

related to the fluid content and fractional volume of pore spaces in subsurface 
layers by expressions such as the following derived from the Patchett-Herrick water 
saturation model for shaly sand formations: ##EQU4## 

Detailed Description Text (24) : 

gives the horizontal resistivity in the reservoir-rock (non-shale) layers of the 
formation. F*.sub.sd,h in equation (5) represents the formation resistivity factor 
for the horizontal resistivity, and B .multidot . Q. sub . v is a factor related to the 
resistivity of "dispersed" shale (shale located within the pore spaces of the 
reservoir rock). V.sub.sh represents the fractional volume within the earth 
formation of interest of the layers of shale ("laminated shale volume"). The other 
terms represent the same quantities as described in the Background section herein. 
See for example, J. G. Patchett et al, "Introduction Section III. Model 
Evaluation", SPWLA Shaly Sand Reprint Volume, Society of Professional Well Log 
Analysts, Houston, Tex. (1982) and M. H. Waxman et al, "Electrical Conductivities 
in Oil Bearing Shaly Sands", SPE Journal, vol. 8, no. 2, pp. 107-122, Society of 
Petroleum Engineers, Richardson Tex., (1968). The porosity can be determined by any 
one of a number of well known measurements, such as acoustic travel time, neutron 
porosity, bulk density, or combinations of measurements such as these as is well 
known in the art. The porosity measurements just described are meant only as 
examples of measurements used to determine the porosity and are not meant to limit 
the invention in any way. 

Detailed Description Text (25) : 

Similarly, for the vertical re s i s t i vi t y z Rt.sub.ver, an expression relating the 
Y3^A£^.L.J3^A^%:^yA^. to the water saturation S.sub.w is: ##EQU5## 

Detailed Description Text (27) : 

In one embodiment of the invention, it has been determined that whatever the values 
of resistivity^ directions, the water saturation, 

S.sub.w, must be equal whether determined from the vertical ^x^hp^jc^^gnt^al 
resistivity measurements. A set of relationships can be developed, where a 
simplifying variable A.sub.w can be defined as: ##EQU6## 

Detailed Description Text (29) : 

where . lambda .. sub . sh represents an "anisotropy factor" relating the y^erjii^a^^and 
horizontal conductivities (or resistivities ) in the same formation. 
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Detailed Description Text (32) : 

Turning now to FIG. 5, a flow chart of the subprocess for determination of water 
saturation according to one embodiment of the invention using the Patchett-Herrick 
model is depicted. Using a logging tool such as disclosed in the '436 patent 
measurements are made within a borehole. The logging tool makes measurements of 
induction signals along and perpendicular to the axis of the instrument as well as 
cross-component signals. As described above, these measurements are processed to 
give a " horizontal "and "vertical" resistivity 220 at each depth in the borehole. A 
measurement of the connate water resistivity R.sub.w is obtained 240. Using the 
value of R.sub.w and a measured value 230 of porosity .PHI. and V.sub.sh, the 
quantities B and Q.sub.v are determined 260 using a relationship given by Juhasz. 
These, together with the values of .PHI. enable the determination of F.sub.horz and 
F. sub. vert using the Waxman-Smit and Patcher-Herrick equations. These are 
substituted in equation (9) along with the measured values of R.sub.tvert and 
R.sub.thorz and an assumed anisotropy factor . lambda .. sub . sh 210 and solved to give 
the value of A.sub.w 250. By substituting this value of A.sub.w in equation (5) 
along with values for B, Q.sub.v and R.sub.w, the water saturation S.sub.w is 
obtained 270. 

US Reference Patent Number (13) : 

6255819 

US Reference Group (13) : 

6255819 20010700 Day et al . 324/303 

CLAIMS : 

I. A method of determining a parameter of interest of a subsurface formation 
containing a sand and a shale surrounding a borehole, the method comprising: (a) 
conveying an electromagnetic logging tool into the borehole and using at least one 
transmitter and one receiver on the tool to obtain measurements indicative of a 
horizontal :„.and .vertical of the formation; (b) obtaining a value for a 
total porosity of said formation, a fractional volume of the shale in said 
formation, and a resistivity of the shale in the formation, (c) using a processor 
for processing said measurements for obtaining a vertical and horizontal 
resistivity of the formation; (d) using a tensor petrophysical model in a processor 
for obtaining from said vertical said obtained total 
porosity, said obtained fractional shale volume and said obtained shale 
resistivity, the parameter of interest wherein said tensor petrophysical model 
comprises a laminar isotropic shale and a laminar sand that is at least one of (i) 

a clean sand, and, (ii) a sand having dispersed shale therein. 

7. The method of claim 1 wherein obtaining a vertical ^ of 
the formation further comprises using the tensor petrophysical model for obtaining 
a laminar shale volume, a laminar sand volume, and a laminar sand conductivity. 

II. The method of claim 1 further comprising using a sensor on the logging tool for 
obtaining an orientation of the tool relative to the subsurface formations and 
using the determined orientation in obtaining said horizontal resistivity and said 
vertical resistivity . 

12. The method of claim 1 further comprises performing an inversion of the obtained 
vertical ajid horizon of the formation . 

13. The method of claim 12 wherein performing an inversion further comprises: (i) 
defining a model including a plurality of layers, each of said plurality of layers 
having a thickness, a horizontals^ a nd a vertical resistivity; and (ii) 
iteratively updating the model based on a difference between an output of the model 
and said obtained horizontal and vertical resistivities. 
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19. A method for determining hydrocarbon saturation of an earth formation 
comprising layers of shale interleaved with layers of reservoir rock, the method 
comprising: (a) determining a vertical resistivity of said earth formation; (b) 
determining a horizontal resistivity of said earth formation; (c) determining a 
porosity of said layers of reservoir rock; and (d) determining a hydrocarbon 
saturation independent of resistivi ty of said layers of shale by combining said 
^^iMi.^sist^it^^ said porosity and said J>o^ 

22. The method as defined in claim 19 wherein said steps of measuring said vertical 
and horizontal resistivity comprise: (a) measuring induction signals along an axis 
of a wellbore drilled through said earth formations; (b) measuring induction 
signals perpendicular to said wellbore axis; (c) measuring cross-component 
induction signals; and (d) determining said vertical and said horizontal 
resistivities using said induction signals measured along and perpendicular to said 
axis and said cross-component induction signals. 

23. A method of obtaining a porosity and water saturation of at least one sand 
layer in a reservoir comprising laminated sands and shales in a subsurface 
formation surrounding a borehole, the method comprising: (a) conveying an 
electromagnetic logging tool into the borehole for obtaining measurements 
indicative of a horizontal and vertical resistivity of the formation; (b) obtaining 
a value for a total porosity of said formation, a fractional volume of the shale in 
said formation, and a resistivity of the shale in the formation; (c) using a 
processor for processing said measurements and obtaining therefrom a layered model 
of vertical and horizontal resistivities of the formation; (d) using a 
petrophysical model for determining from the layered model of vertical and 
horizontal ..resistiviti the porosity and water saturation of said sand. 

24. The method of claim 23 wherein obtaining a resistivity of shale further 
comprises conveying a galvanic tool having button electrodes thereon into the 
borehole and averaging in a vertical direction and in a horizontal direction 
measurements made by said button electrodes. 

28. A method of determining a parameter of interest of a subsurface formation 
containing a sand and a shale surrounding a borehole, said subsurface formation 
having a total porosity, a fractional volume of the shale in said formation, a 
resistivity of the shale, a vertical .re the 
method comprising using a tensor petrophysical model in a processor for obtaining 
from said vertical and horizontal resistivities, said obtained total porosity, said 
obtained fractional shale volume and said obtained shale resistivity, the parameter 
of interest wherein said tensor petrophysical model comprises a laminar isotropic 
shale and a laminar sand that is at least one of (i) a clean sand, and, (ii) a sand 
having dispersed shale therein. 
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L15: Entry 3 of 3 



File: USPT 



Jul 3, 2001 



DOCUMENT- IDENTIFIER : US £255819 Bl 

TITLE: System and method for geologically-enhanced magnetic resonance imaging logs 



Brief Summary Text (9) : 

For example, many of the low resistivity and low contrast pay sands are known to be 
thinly laminated sand/silt or shale sequences. Most commercially available 
resistivity logging tools have coarse vertical resolution and fail to read 
resistivity of individual sand or shale layers. Instead, they read only the 
averaged horizontal resistivity that is low and dominated by the high conductivity 
of silt and shale layers, although individual sand layers can be highly resistive. 
As a result, the low resistivity of the layer sequence may be incorrectly 
interpreted as poor hydrocarbon potential of the formation. 
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PGPUB- DOCUMENT -NUMBER : 20040010373 
PGPUB-FILING-TYPE : new 

DOCUMENT-IDENTIFIER: US 20040010373 Al 

TITLE: INTELLIGENT DIAGNOSIS OF ENVIRONMENTAL INFLUENCE ON WELL LOGS WITH MODEL- 
BASED INVERSION 

PUBLICATION-DATE: January 15, 2004 



INVENTOR- INFORMATION: 








NAME 


CITY 


STATE 


COUNTRY 
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□ 2. Document ID: US 6005389 A 

L29: Entry 2 of 2 File: USPT Dec 21, 1999 

US-PAT-NO: 6005389 

DOCUMENT-IDENTIFIER: US 6005389 A 

TITLE: Pulse sequences and interpretation techniques for NMR measurements 
DATE-ISSUED: December 21, 1999 
INVENTOR- IN FORMAT I ON : 

NAME CITY STATE ZIP CODE COUNTRY 

Prammer; Manfred G. Downingtown PA 

US-CL-CURRENT: 324/303; 324 / 300 
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